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Abstract 

We derive in-medium nucleon-nucleon (NN) cross sections in a microscopic 
model. Our calculations are based upon the Bonn NN potential and the 
Dirac-Brueckner approach for nuclear matter. We consider energies up to 
300 MeV (in the laboratory frame) and densities up to twice nuclear matter 
density. Our results deviate substantially from cross section parametrizations 
that are commonly used in the nuclear medium. 
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I. INTRODUCTION 



An exciting topic in contemporary nuclear physics is the study of the medium (density 
and/or temperature) dependence of the properties of hadrons and hadronic processes. A 
well-known example is the mass of the nucleon which decreases with increasing density. This 
is implied by relativistic (e.g. Walecka model and Dirac-Brueckner approach P-pi]|]) as 
well as non-relativistic (e.g. momentum dependent Skyrme forces |[TT|JT^) models. Based on 
effective chiral Lagrangians, Brown and Rho [0,0 argue that the masses of other hadrons 
also decrease with increasing density, at about the same rate as the mass of the nucleon. 
Recent finite density QCD sum rule calculations also indicate that the masses of hadrons 
should decrease with increasing density |T5Hr7|. 

Not only the static properties of hadrons (e.g. their masses) but also the dynamical ones 
(e.g. two-body scattering ) are modified due to the presence of the medium. In-medium 
nucleon-nucleon (NN) scattering differs from the corresponding one in free space, mainly due 
to Pauli blocking of intermediate and final states as well as the mean field. In conjunction 



with nucleus-nucleus collisions, Faessler and co-workers |18|-p0| have studied in-medium NN 
scattering, based on non-relativistic Brueckner calculations and the Reid soft-core potential. 
This approach has been applied in calculations of the nucleus-nucleus optical potential at 
low energies [|T^-^ and in the transport model description of nucleus-nucleus collisions 



at intermediate energies ||24-P4|. It should, however, be emphasized that, non-relativistic 
many-body theory, e.g. the Brueckner approach p5H27|| and the variational method ||28H30|, 
is not able to reproduce correctly the saturation properties (saturation density and energy) 
of nuclear matter, when two-body forces are applied 

The investigation of in-medium NN scattering is of interest for intermediate-energy 
heavy-ion reactions. Experimentally, nucleus-nucleus collisions at intermediate energies pro- 
vide a unique opportunity to form a piece of nuclear matter in the laboratory with a density 
up to 2-3po (with po, in the range of 0.15 to 0.19 fm~'^, the saturation density of normal 



nuclear matter; in this paper we use po=0.18 fm ^) [pT],p2|. Thus it is possible to study 
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the properties of hadrons in dense medium. Since this piece of dense nuclear matter exists 
only for a very short time (typically 10~^^ -10^^^ s), it is necessary to use transport models 
to simulate the entire collision process and to deduce the properties of the intermediate 
stage from the known initial conditions and the final-state observables. At intermediate 
energies, both the mean field and the two-body collisions play an equally important role in 
the dynamical evolution of the colliding system; they have to be taken into account in the 
transport models on an equal footing, together with a proper treatment of the Pauli blocking 
for the in-medium two-body collisions. The Boltzmann-Uehling-Uhlenbeck (BUU) equation 
||33| , |3^ and quantum molecular dynamics (QMD) |3^j3^, as well as their relativistic exten- 



sions (RBUU and RQMD) [|T^,^^-^, are promising transport models for the description 
of intermediate-energy heavy-ion reactions. In addition to the mean field, the in-medium 
NN cross section is also an important ingredient of these transport models. Specifically, 
in-medium total as well as differential NN cross sections are needed by these models in 
dealing with the in-medium NN scattering using a Monte Carlo method. Moreover, if these 
models are used to calculate the spectrum of particles produced in heavy-ion collisions, such 
as photons, pions and kaons, one also needs to know the in-medium elementary particle 
production cross section e.g., the in- medium cross section of neutron-proton 

(np) bremsstrahlung {np np'y) is needed in the calculation of the photon production cross 
section in nucleus-nucleus collisions. Up till now, most calculations of particle production 
in nucleus-nucleus collisions are based on the free production cross sections. 

It is the purpose of this paper to calculate the elastic in-medium NN cross sections in a 
microscopic way. We base our investigation on the Bonn meson-exchange model for the NN 
interaction |]5|,^ and the Dirac-Brueckner approach 0-0 for nuclear matter. 



Speaking in more general terms, it is the fundamental goal of traditional nuclear physics 
to describe nuclear structure and nuclear reactions in terms of the same realistic NN inter- 
action. This bare interaction should be based as much as possible on theory and describe 
the two-nucleon system (NN scattering as well as the deuteron) accurately. 

There are two aspects to the problem. First, one needs a realistic NN interaction which 



is ultimately determined by the underlying dynamics of quarks and gluons and should in 
principle be derived from quantum chromodynamics (QCD). However, due to the nonper- 
turbative character of QCD in the low-energy regime relevant for nuclear physics, we are 
far away from a quantitative understanding of the NN interaction in this way. On the other 
hand, there is a good chance that conventional hadrons, like nucleons and mesons, remain 
the relevant degrees of freedom for a wide range of low-energy nuclear physics phenomena. 
In that case, the overwhelming part of the NN potential can be constructed in terms of 
meson-baryon interactions. In fact, the only quantitative NN interactions available up till 
now are based upon the idea of meson-exchange; two well known examples are the Paris 
potential ||46|J47[] and the Bonn potential |]5|,^ (see Ref. for a comprehensive overview of 
the history of meson-exchange models). 

The second aspect of the problem concerns a suitable many-body theory that is able to 
deal with the bare NN interaction which has a strong repulsive core. The Brueckner approach 
25|^7[ and variational method p8|-p0| have been developed for this purpose. However, when 



using two-body forces, the Brueckner approach and the variational method are not able to 
reproduce correctly the saturation properties of nuclear matter. Inspired by the success 
of the Dirac phenomenology in intermediate-energy proton-nucleus scattering p5|-pD|] and 



the Walecka model (QHD) for dense nuclear matter |[T|J^, a relativistic extension of the 
Brueckner approach has been initiated by Shakin and co-workers 0,^, frequently called 
the Dirac-Brueckner approach. This approach has been further developed by Brockmann 
and Machleidt P,|6|,|5T| and by ter Haar and Malfliet pHlOO. Formal aspects involved in the 
derivation of the relativistic G-matrix have been discussed in detail by Horowitz and Serot 



| 5^ , |53[| . The common feature of all Dirac-Brueckner results is that a repulsive relativistic 
many-body effect is obtained which is strongly density-dependent such that the empirical 
nuclear matter saturation can be explained. The Dirac-Brueckner approach thus provides a 
reasonable starting point for pursuing the longstanding goal of self-consistently describing 
nuclear matter, finite nuclei and nuclear reactions based on the same realistic NN interaction. 
The extension of this relativistic approach to the domain of finite nuclei has been discussed 



by Miither et al. |Q and recently by Brockmann and Toki |5^, while the investigation of 
nucleon-nucleus scattering has been initiated by MacKellar et al. [Q. 

This paper is the first in a series in which we investigate systematically in-medium NN 
scattering based on the Bonn potential and the Dirac-Brueckner approach. In this paper, we 
are concerned with elastic in-medium NN scattering which is the most important two-body 
process in nucleus-nucleus collisions at incident energy below 300 MeV per nucleon. In Sec- 
tion 2 we give a brief description of the Bonn potential and compare theoretical predictions 
with experimental data for free-space NN scattering. The Dirac-Brueckner approach and 
the predictions for nuclear matter are discussed in Section 3. The results for the in-medium 
NN cross sections are presented and discussed in Section 4. Finally we give a brief summary 
and outlook in Section 5. 

II. THE BONN MODEL AND NN OBSERVABLES 



Two-nucleon scattering is described covariantly by the Bethe-Salpeter equation . As 
this four- dimensional integral equation is very difficult to solve, so-called three-dimensional 
reductions have been proposed, which are more amenable to numerical solution ||58| , |59| . 
Among the different forms of three-dimensional reductions, the one suggested by Thompson 
is particularly suitable for the relativistic many-body problem 1^,^. In terms of the 



i?-matrix (or fT-matrix) the Thompson equation reads in the center-of-mass system ||6^ 



<,) = VW, q) + P / I^J^q'. q) (1) 

where q, k and q' are initial, intermediate and final relative momenta, resp., of the two 
scattering nucleons. Eq = {m? + q^)^/^ with m the mass of the free nucleon. V denotes the 
principal value. 

In the one-boson-exchange (OBE) model, the kernel of this integral equation, V^(q', q), is 
the sum of one-particle-exchange amplitudes of certain bosons with given mass and coupling. 
In the OBE Bonn model P, six nonstrange bosons with mass below 1 GeV are used; they 
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are pseudoscalar mesons (vr and 77), scalar mesons {S and a) and vector mesons (p and u). 
The meson-nucleon interactions are described by the following Lagrangians 

Cp, = -Il^^^^Y^pd^cpiP^^ (2) 
Cs = gsi'n^''^ (3) 



with il) the nucleon and ^["n the meson fields (notation and convention as in Ref. ||6T]]). For 



isospin-1 mesons (vr, 5 and p ), is to be replaced by r ■ with r the usual Pauli 
matrices. Note that the pseudovector coupling is used for pseudoscalar mesons in order to 
avoid unrealistically large antiparticle contributions. 

From the above Lagrangians, we can derive the OBE amplitudes; e.g., the contribution 
from the isoscalar-scalar meson is given by 

< q'A;A^|V;«^^|qAiA2 >= -glu{c(, A;)«(q, Ai)u(-q', A^)«(-q, ^2)^-^^^^^^^ (5) 

\(\ qj + nig 

where Ai (A^), with i = 1,2, denotes the helicity of the incoming (outgoing) nucleons. 
jFs[(q' — q)^] is a form factor of monopole type which simulates the short-range physics 
governed by quark-gluon dynamics: 

Aj + (q' - q)^ 

with As the cutoff mass of the isoscalar-scalar meson. The Dirac spinors are normalized 
covariantly 

u{q,X)u{q,X) = l (7) 



The OBE amplitudes for other mesons are given in Refs. P J^ , |60[] . 

Three sets of OBE potential parameters, denoted by Bonn A, B and C, have been 
proposed We reprint them in Table 1. The main difference between the three parameter 



sets is the cut-off parameter for the irNN vertex, which is 1.05, 1.2 and 1.3 GeV for Bonn 
A, B and C, respectively. Consequently, the three potentials differ in their strength of the 
tensor force component which is not well constrained by present NN data. Bonn A has the 
weakest tensor force and Bonn C the strongest. All three potentials reproduce the deuteron 
properties and the phase shifts of NN scattering up to about 300 MeV accurately (cf. Refs. 



Our goal is to calculate the in-medium NN cross sections microscopically. For this 
purpose it is important that the bare NN interaction describes the free NN scattering cross 
sections correctly. In Fig. 1, we compare the predictions by the Bonn potentials with the 
neutron-proton {np) differential cross section data at 50, 129 and 212 MeV. It is seen that all 
three potentials reproduce the data very well. In Fig. 2 we show the np total cross sections 
in the energy range 50-300 MeV. There is good agreement between theory and experiment. 

III. DIRAC-BRUECKNER APPROACH AND NUCLEAR MATTER 

PROPERTIES 

As mentioned in the Introduction, the essential point of the Dirac-Brueckner approach is 
the use of the Dirac equation for the description of the single-particle motion in the nuclear 
medium 



where Us is an attractive scalar field and Uy the timelike component of a repulsive vector 
field, m is the experimental mass of the free nucleon. The solution of this equation is 



IS- 



[q: • k P{m + Us) + Uv]u{k, s) = eku{k, s) 



(8) 





with 



m = m + Us 



Efc = (m2^k2)V2 
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and Xs is a Pauli spinor. Note that the in-medium Dirac spinor is obtained from the free 
Dirac spinor by replacing m by m. 

The single-particle energy resulting from Eq. (8) is 

ek = Ek + Uv. (10) 

Similar to conventional Brueckner theory, the basic quantity in the Dirac-Brueckner 
approach is a G-matrix which satisfies the in-medium Thompson equation (also known as 
relativistic Bethe-Goldstone equation) 1^-0, which reads in the nuclear matter rest frame: 

G(q',q|P,5)=^(q',q)+P / -f^^ (q', k) ^;^^^^fc^G(k, q|P, 5) (11) 

l^^J -^(l/2)P+k ^ ^ ^-^(l/2)P+k 

with 

Z = 2£^(i/2)p+q 

and P is the cm. momentum of the two colliding nucleons in the nuclear medium. Equa- 
tion (11) is density-dependent which is suppressed in our notation. Notice that in the 
energy denominator of Eq. (11) (which is the difference of single-particle energies of the 
kind Eq. (10)), Uy drops out since it is constant. 

The in-medium Thompson equation differs from the free one mainly in three points 
(comparing eq. (11) with eq. (1)). First, the Pauli operator Q prevents scattering into 
occupied intermediate states ('Pauli effect'). Note that this is different from the Pauli 
blocking factor for the final states which is always included in the transport models describing 
nucleus-nucleus collisions. Second, the nucleon mean field due to the medium reduces the 
mass of the nucleon and affects the energy denominator in eq. (11) which is now density 
dependent, while in eq. (1) the energy denominator uses free relativistic energies ('dispersion 
effect'). Finally and most importantly, the potential used in the in-medium Thompson 
equation, as indicated by the tilde, is evaluated by using the in-medium Dirac spinors of eq. 
(9) instead of the free ones that are used for the \^ in eq. (1). This leads to the suppression 
of the attractive a exchange which increases with density. The fact that the Dirac-Brueckner 
approach is able to reproduce quantitatively the saturation properties of nuclear matter is 
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mainly due to this relativistic effect. This observation also implies that the in-medium NN 
cross sections based on the non-relativistic Brueckner approach lack one important aspect 
20| , |22|p3| , namely, the effect which is due to the medium modification of the potential. 



The scalar and vector fields of the Dirac Eq. (8) are determined from 



in . ^ im? 



Us + Uy= Y: ^<^J\GmJ-J^> (12) 

-C/j j<kF ^iJ^j 

which is the relativistic analogue to the non-relativistic Brueckner-Hartree-Fock definition 
of a single-particle potential. 

Since the kernel of the in-medium Thompson equation (eq. (11)) depends on the solution 
of the Dirac equation (eq. (8)), while for the Dirac equation one needs the scalar and vector 
potentials which are related to the G matrix via Eq. (12), one has to carry out an iterative 
procedure with the goal to achieve self-consistency of the two equations 0-0 : starting from 
reasonable initial values for uf^ and Uy^ , one solves the in- medium Thompson equation in 



momentum space by means of the matrix inversion method |7T| to get the G matrix which 
leads by means of Eq. (12) to a new set of values for U^^ and Uy^ to be used in the next 
iteration; this procedure is continued until convergence is achieved. 

The nuclear equation of state, that is the energy per nucleon, S/A, as a function of 
density, p, is obtained from the G-matrix P,^] : 

S 1 ^ mrh + p'^ 1 . ~ ,, , . 

A = A^ F ^ 2A ^ jjr < ^J\G{z)\^J - J^ > -m (13) 

^ ^ i<kp i,j<kF ^i^j 

We show in Fig. 3 the energy per nucleon S/A as function of nuclear matter density, 
p/Po (Po=0.18 fm~^). The solid, dashed and dotted lines are the results corresponding to the 
Bonn A, B and C potentials, respectively. The open rectangle indicates the empirical region 
for nuclear matter saturation. There is some difference between the results obtained with the 
different potentials, especially for the nuclear equation of state at higher densities, although 
these potentials predict the same in free-space NN scattering. This difference can be traced 
back to differences in the strength of the tensor force, as reflected in the predictions for the 
D-state probability, Pd, of the deuteron |]5|, |45| , |70| . The Bonn A potential, with the lowest 
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D-state probability (Pd=4.47%), indicating a weak tensor force, yields the best prediction 
for the empirical values of nuclear matter saturation. More results and discussion of the 
properties of (dense) nuclear matter and neutron matter can be found in Refs. 



IV. IN-MEDIUM NN CROSS SECTIONS 

In the previous sections, we discussed briefly the Bonn potentials, their predictions for 
free-space NN scattering and the Dirac-Brueckner approach for nuclear matter. The poten- 
tials describe free-space NN scattering accurately and nuclear matter saturation is repro- 
duced well in the Dirac-Brueckner approach (especially with the Bonn A potential). 

In-medium NN cross sections can be calculated directly from the G matrix p2|,|9|,|69 



Alternatively, one can also calculate first the in-medium phase shifts, which are defined in 
terms of the partial-wave G-matrix elements like the free-space NN phase shifts are defined 



in terms of the i?-matrix elements. From the in-medium phase shifts, the in-medium 
NN cross sections are obtained in the usual way. 

We calculate the G-matrix, from which we obtain the in-medium cross sections, in the 
center-of-mass (cm.) frame of the two interacting nucleons, i. e., we use Eq. (11) with 
P = 0. For the starting energy in Eq. (11), we have now: z = 2Eq = 2y/mP~+~q^, where 
q is related to the kinetic energy of the incident nucleon in the "laboratory system" (Eiab), 
in which the other nucleon is at rest, by: Eiab = 2q^/m. Thus, we consider two colliding 
nucleons in nuclear matter. The Pauli projector is represented by one Fermi sphere as in 
conventional nuclear matter calculations. This Pauli-projector, which is originally defined 
in the nuclear matter rest frame, must be boosted to the cm. frame of the two interacting 
nucleons. For a detailed discussion of this and the explicit formulae, see Refs. [|8|,|53|. In 
summary, for in- medium NN scattering i^'-matrix, we use the G-matrix of Eq. (11) with 
P = 0, 5 = 2Eq, and the Pauli-projector Q replaced by the ellipsoidal one due to Lorentz 
boosting. 

In Fig. 4 we show the in-medium np phase shifts for the ^S'o, ^Pq and "^Si partial-wave 
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state as a function of the laboratory energy. For each partial waves, we present three different 
results corresponding to density p=0 (free-space scattering, solid lines), po (dashed lines), 
and 2po (dotted lines). The results presented in this figure are obtained with the Bonn A 
potential. A clear decrease of the NN phase shifts with increasing density is observed. This 
is due to the Pauli and dispersion effect as well as the relativistic medium effects. 

In Fig. 5 we show the results for the in-medium np differential cross sections as a function 
of the cm. angle for Eiab—50, 100, 250 and 300 MeV. For each incident energy, we present 
three different results corresponding to the medium density p—0 (free-space scattering, solid 
lines), (l/2)po (long-dashed lines), po (short-dashed lines) and 2po (dotted lines). The results 
are obtained by using the Bonn A potential. At low incident energies (Figs. 5(a) and 5(b)), 
the np differential cross section always decreases with increasing density, at both forward 
and backward angles. At high incident energies (Figs. 5(c) and 5(d)), the np differential 
cross section for forward angles decreases when going from p—0 to (l/2)po and then increases 
for higher densities. The differential cross section at backward angles always decreases with 
density. While the free np differential cross sections are highly anisotropic, the in-medium 
cross sections become more isotropic with increasing density. 

In Fig. 6, we compare the in-medium np differential cross sections as predicted by 
different potentials. The cross sections in free space are shown in Fig. 6(a), while Fig. 
6(b) and 6(c) display them for p — pq and p — 2po, respectively. The incident energy is 
fixed at 100 MeV. While the predictions for the np differential cross section in free space 
are essentially the same, these potentials lead to slightly different in-medium np differential 
cross sections at higher densities (see Fig. 6(b) and 6(c)): at p — po, the prediction by the 
Bonn A potential is larger than those by Bonn B and C in all directions, while at p = 2po, 
the prediction by the Bonn A potential is smallest in the forward direction and largest in 
the backward direction. This difference is mainly due to differences in the strength of the 
tensor potential. The purpose of Fig. 6 is to give some idea of the model dependence in our 
approach, due to our incomplete knowledge of the nuclear force. Fortunately, this model 
dependence turns out to be quite moderate. 
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It is also interesting to compare the results obtained in this work with the parametrized 



NN cross sections proposed by Cugnon et al. p3| , [72| which are often used in transport mod- 



els such as BUU and QMD. This is done in Fig. 7 for a laboratory energy of 300 MeV. 
We compare our predictions for three different densities (using the Bonn A potential) with 
the Cugnon parametrization. It is clearly seen that, while the Cugnon parametrization is 
almost isotropic, the microscopic results still have some anisotropy at all densities consid- 
ered. The anisotropy in the present results decreases with increasing density (mainly due 
to the decrease of the magnitude of the ^Pi phase shift). There is also density dependence 
in the microscopic differential cross section, while the Cugnon parametrization is density 
independent. 



We mention that Cugnon et al. |72| have parametrized the free pp cross sections. This 



explains the almost-isotropy in their differential cross sections as well as the lack of the 



density dependence. Note that in the work of Cugnon et al., [^], no difference is made 
between proton and neutron; thus, the pp cross sections are also used for np scattering. 
There are, however, well-known differences between pp and np cross sections which in the 
more accurate microscopic calculations of the near future may be relavant. The difference 
between the (in-medium) np and pp cross sections will be discussed in detail in a forthcoming 
paper ||73|. In this paper, all "NN cross sections" are np cross sections. 



The np differential cross section in free space can be well parametrized by the following 
simple expression: 

^{Eiah, e) = ^^4?7 -expMcos'd + sin"- - 1.0)] (14) 

dn^ ' ' 1.0 + 0.05(Efi - 15.5) ^ 1 ^ ' 



with 



and 



0.0008(E°;f - 4.625) for Eiab < lOOMeV 



b = 0.0006(36.65 - ) for Eiab > lOOMeV. 
with Eiab in the units of MeV. 
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The quality of this parametrization is demonstrated in Fig. 8, where it is compared with 
the microscopic results (points) based on Bonn A at three energies. It would be useful to 
parametrize the in-medium np differential cross section as well. However, the complicated 
dependence of the in-medium differential cross sections on angles, energy, and especially 
density makes this very difficult. Instead, we have prepared a data file, containing in-medium 
differential cross section as a function of angle for a number of densities and energies, from 
which the differential cross sections for all densities in the range 0-3po and all energies in 
the range 0-300 MeV can be interpolated. This data file is available from the authors upon 
request. 

In addition to the in-medium NN differential cross sections which enter the transport 
models to determine the direction of the outgoing nucleons, the in-medium total NN cross 
sections, ct^n, are also of interest. They provide a criterion for whether a pair of nucleons 



will collide or not by comparing their closest distance to \J(Jnn/t^- We show in Fig. 9(a) 
the in- medium total cross sections as function of the incident energy Eiab and in Fig. 9(b) 
as function of density. For completeness, we also list in Table 2 the microscopic in-medium 
np total cross sections as function of density and energy. All results are obtained by using 
the Bonn A potential. 

It is seen that the in-medium total cross sections decrease substantially with increasing 
density, particularly for low incident energies. This agrees well with the findings of Ref . , 
but disagrees with Ref p2[. In Ref. [^, an enhancement of the cross section was found 



below 150 MeV and little change (as compared to the free cross section) above 150 MeV. The 
major difference between our (and Ref 's. ||^) calculations and the ones of Ref. p2[ is that we 



include relativistic many-body effects, which are ignored in Ref. p2|]. Besides this, we as well 



as Ref. include the Pauli and dispersion effects (which are larger than the relativistic 
effects). The latter two effects reduce considerably the magnitude of, in particular, the 
^5*0 and ^Si G-matrix elements. This fact is well-known since the work of Bethe and his 
group on Brueckner theory in the 1960's WM- This leads to a substantial reduction of the 



binding energy in nuclear matter as well as the in-medium NN cross sections, since both are 
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based upon the G-matrix. On the background of these well-estabhshed facts, it is hard to 



understand the results of Ref. [Q, while our results and the ones of Ref. are physically 
quite meaningful. 

At higher energies (~ 300 MeV), the medium effect becomes smaller as compared to 
lower energies (cf. Fig. 9a), but it does not vanish. We use the continuous choice for the 
single-particle potential, and so there is a dispersive effect also for higher momenta; in 
addition, there are the relativistic effects. This explains the non-negligible medium effects 



at intermediate energies. Again, in this we agree with Ref. 0, and disagree with Ref. ||22 
where above 150 MeV the free cross sections were obtained. The use of the 'gap-choice' for 
the single-particle spectrum (i. e., free energies above the Fermi surface) and the omission 
of relativistic effects in Ref. p2| may be the explanation here. 

In Fig. 10 we compare the total in-medium cross sections obtained with different poten- 
tials. Figs. 10(a) and 10(b) correspond to the density p = po and 2po, respectively. The 
solid, dashed and dotted lines refer to Bonn A, B and C, respectively. By comparing Fig. 10 
with Fig. 2 we see that, although the potentials lead to essentially the same predictions in 
free space, there is some difference in the medium, and this difference increases with density. 
The prediction by the Bonn A potential is the largest at low energies and the smallest at 
high energies. Again, forturnately, the model dependence is moderate. 

We compare in Fig. 11 the total in- medium cross section obtained in the present work 



using the Bonn A potential with the one used by Cugnon et al. [|33| , |72| . The solid, long- 
dashed and short-dashed lines are the present results corresponding to the medium den- 
sity p = 0, l/2po and 3/2po, respectively, while the dotted line represents the Cugnon 
parametrization. It is seen that at low energies and low densities, the Cugnon parametriza- 
tion underestimates the microscopic results, while at higher energies and higher density it 
is the other way around. Note that the Cugnon parametrization is not density dependent, 
and, thus, predicts the same for all densities. 

Finally, we propose a parametrization for the total np cross section as a function of the 
incident energy Eiab and density p 
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arUEia,, p) = (31.5 + 0.092abs(20.2 - ^L'')"^)^^^;^^^^ (15) 

where Ei^b and p are in the units of MeV and fm~^, respectively. 

We compare in Fig. 12 this parametrization (hnes) with the microscopic resuhs (points) 
based on Bonn A at four densities. It is seen that the parametrization reproduces well the 
microscopic results for energies and densities considered in this work. 



V. SUMMARY AND OUTLOOK 

In this paper, we presented a microscopic derivation of elastic in-medium NN scattering 
cross sections for energies up to 300 MeV and densities up to 2po- This investigation is based 
upon the Bonn NN potential and the Dirac-Brueckner approach for nuclear matter. 

The 'bare' Bonn potential reproduces the free NN scattering cross sections (differential 
as well as total) accurately. When the potential is applied to nuclear matter using the 
Dirac-Brueckner approach, the saturation properties (saturation density and binding energy) 
are reproduced correctly. Thus, the Bonn model provides a good starting point for other 
investigations. The major conclusions of the present microscopic calculations are: 

(1) There is strong density dependence for the in-medium cross sections. With the 
increase of density, the cross sections decrease. This indicates that a proper treatment of 
the density-dependence of the in- medium NN cross sections is important. 

(2) Our microscopic predictions differ from the commonly used parametrizations of the 
differential and the total cross sections developed by Cugnon et al. p3| , [72| . The Cugnon 



parametrization underestimates the anisotropy of the in-medium np differential cross sec- 
tions. In the case of the total cross sections, the Cugnon parametrization either underesti- 
mates or overestimates the microscopic results, depending on energy and density. 

At energies above 300 MeV, inelastic channels enter the picture. Microscopic models 
that also describe the inelasticity have to be applied. These models will then also allow 
to calculate the in-medium pion production. This is under investigation. 
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Table 1 

Parameters of the OBE Bonn potentials used in this work (reprinted from Ref. 
Appendix A, Table A. 2, therein). 

Bonn A Bonn B Bonn C 



mesons m„ (MeV) gH^n K (GeV) gl/Ari K (GeV) gl/Ar: K (GeV) 





138.03 


14.9 


1.05 


14.6 


1.2 


14.6 


1.3 




548.8 


7 


1.5 


5 


1.5 


3 


1.5 




769 


0.99 


1.3 


0.95 


1.3 


0.95 


1.3 




782.6 


20 


1.5 


20 


1.5 


20 


1.5 


6 


983 


0.7709 


2.0 


3.1155 


1.5 


5.0742 


1.5 


a 


550 


8.3141 


2.0 


8.0769 


2.0 


8.0279 


1.8 


a) ga = 


tjc^ for 


vr and t]. 













We use ^=6.1 and ^=0. 

9p 
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Table 2 

Microscopic in-medium np cross section atot (mb) obtained in this work with Bonn A. 
Po=0.18 fm-^ 



Eiab (MeV) 



p 


50 


100 


150 


200 


250 


300 





164.8 


72.18 


49.17 


39.57 


34.27 


30.85 


(l/2)po 


93.92 


47.96 


26.71 


20.62 


18.76 


18.14 


Po 


54.67 


28.50 


18.17 


15.14 


14.93 


15.40 


(3/2)po 


37.76 


20.93 


15.53 


13.91 


14.41 


15.26 


2po 


28.48 


17.96 


15.96 


14.94 


15.14 


15.97 
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FIGURES 

FIG. 1. Neutron-proton differential cross sections at (a) 50 MeV, (b) 129 MeV and (c) 212 
MeV. The curves are predictions by the Bonn potentials. The data at 50 MeV arc from Ref. [62] 
(solid circles) and Ref. [63] (solid squares), while the data at 129 MeV and 212 MeV are from Ref. 
[64] and Ref. [65], respectively. 

FIG. 2. Neutron-proton total cross sections in the energy range 50-300 MeV. The curves are 
the predictions by the Bonn potentials. The data are from Ref. [66] (solid circles), Ref. [67] (solid 
triangles) and Ref. [68] (solid squares). 

FIG. 3. Energy per nucleon in nuclear matter as obtained in the Dirac-Brueckner approach 
using the Bonn potentials. The box denotes the empirical region of nuclear matter saturation. 

FIG. 4. In-medium np phase shifts for (a) ^Sq, (b) ^Pq and (c) ^S\. The results are obtained 
with the Bonn A potential. 

FIG. 5. In-medium np differential cross sections at (a) 50 MeV, (b) 100 MeV, (c) 250 MeV and 
(d) 300 MeV laboratory energy, as obtained for various densities. The Bonn A potential is used. 

FIG. 6. In-medium np differential cross sections at 100 MeV laboratory energy for the densities 
(a) p=0, (b) p = po and (c) p = 2po- Predictions by three differential potentials are shown. 

FIG. 7. In-medium np differential cross section at 300 MeV laboratory energy for the densities 
p = 0, p = Po and p = 2po, using the Bonn A potential. The dotted curve is the Cugnon 
parametrization [33,72] 

FIG. 8. Free-space np differential cross sections at three energies. The lines are the predictions 
by the parametrization, eq. (14), the squares, circles and triangles are the microscopic results at 
Eiab=50, 100 and 300 MeV, respectively, as predicted by the Bonn A potential. 

FIG. 9. In-medium np total cross sections (a) as function of the incident energy and (b) as 
function of density. The results are obtained with the Bonn A potential. 
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FIG. 10. In-medium np total cross sections as predicted by different potentials at densities (a) 
p = po and (b) p = 2po. 

FIG. 11. In-medium np total cross sections as described by the Cugnon parametrization [33,72] 
(dotted line) are compared with the predictions by our microscopic calculation using the Bonn A 
potential. 

FIG. 12. In-medium np total cross section as function of incident energy at four densities. The 
curves are the description by the parametrization, eq. (15), the squares, circles, triangles and 
diamonds are the microscopic results at p = 0, (l/2)po, po and 2po, respectively, using the Bonn 
A potential. 
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The figures are available upon request from 

GQLI@TAMCOMP.BITNET 

Please, include your mailing address with your request. 
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